Interface-Localized Mode In Bilayer Film 
Ferromagnetic Resonance Spectrum* 

H. Puszkarski** and S. Mamica 

Surface Physics Division, Faculty of Physics, Adam Mickiewicz University, 
61-614 Poznah, ul. Umultowska 85, Poland 

PACS numbers: 75.70.-i; 76.50.+g; 75.30.Ds 

*To be published in "Horizons in World Physics" Vol. 243 (Nova Science Publishers) 
** Corresponding author: henpusz@main.amu.edu.pl 

ABSTRACT 

Ferromagnetic resonance (FMR) in exchange-coupled bilayer films has been the 
subject of intensive studies in recent years. From the experimental viewpoint, a char- 
acteristic feature of this FMR is that some specimens show single resonance, whereas 
others show double resonance. Moreover, double resonance can exhibit a regular pat- 
tern, in which the high-field (HF) line intensity surpasses that of the low-field (LF) 
line, or it can exhibit an inverted pattern with the HF line less intense than the LF 
line. There is a general agreement that the inverted FMR pattern occurs when the HF 
line is an 'optic mode ', i. e. an out-of-phase composition of individual sublayer modes, 
and the LF line is an 'acoustic mode', or an in-phase mode composition. The exist- 
ing theoretical explanations of bilayer ferromagnetic resonance are, as a rule, based 
on phenomenological equation of motion of the magnetization vector. In this paper, 
we propose a theory of FMR in ultrathin bilayers based on an entirely microscopic 
approach, using the Heisenberg model of localized spins and assuming that the two fer- 
romagnetic sublayers are exchange-coupled through their interface. Both the strength 
and the sign of this interface coupling {Jint) is arbitrary (we admit ferromagnetic or 
antiferromagnetic interface coupling). The Hamiltonian contains an exchange energy 
and a Zeeman energy terms; the external field is assumed to be applied obliquely to the 
film surface. We focus on the effects originating from the interface coupling, though the 
system is assumed to exhibit also pinning effects, originating from surface anisotropy 
on the outer surfaces of the film as well as from intrinsic interface anisotropy present 
on internal interfaces. The latter anisotropy is assumed to consist of two components: 
uni-directional {Kint) and uni-axial [Dint)- We show that the resonance spectrum in 
symmetric bilayer is completely independent of Jmt, but depends strongly on the ap- 
plied static field configuration with respect to the interface normal (angle 9). A critical 
angle Ocrit is found to exist (as in the case of single-layer film) for which the multipeak 
spectrum reduces to a single-peak one. This rigorous microscopic FMR theory does 
explain the inverted pattern of the bilayer FMR spectrum by assuming the HF line to 
correspond to an in-phase mode, hut of interface-localized nature. The intensity of such 
localized mode decreases with growing strength of its localization at the interface and, 
when the localization becomes sufficiently strong, becomes lower than the intensity of 
the other mode (which is of the bulk type). This gives a possibility to explore the HF 
resonance line corresponding to the interface-localized mode as a potential source of 
information concerning the bilayer interface. 
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1 Introduction 

Magnetic multilayers have recently become the subject of intensive studies, both 
theoretical and experimental, in which special attention is paid to interface parameters, 
such as interface exchange coupling or interface anisotropy. This interest in interface 
is due to its significant effect on the properties of the multilayer system as a whole. 
One of the key methods of interface investigation is based on ferromagnetic resonance 
(FMR). The interface coupling as well as the interface anisotropy |31 !SI in a 

multilayer can be studied by means of the FMR spectrum. In anisotropy studies, 
angle relations play an important role (HJ-IEj- Also, FMR has been recently used in 
investigating magnetic particles size distribution il9|. 

From the experimental viewpoint, the FMR in bilayer films is characterized by the 
fact that some specimens show single resonance, whereas others show double resonance. 
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Moreover, the double resonance can exhibit a regular pattern, in which the high-field 
(HF) line possesses an intensity greater than that of the low-field (LF) line, or it can 
exhibit an inverted pattern with a HF line less intense than the LF line. According to 
a commonly accepted interpretation, the inverted FMR pattern occurs when the HF 
line is an out-of-phase composition of the individual sublayer modes {"optic mode") 
and the LF line is an in-phase mode composition {"acoustic mode"). Here, consensus 
seems to be due to the circumstance that out-of-phase modes, naturally, are associated 
with lower net magnetization; this is invoked as an explanation of the lower intensity 
of the HF resonance line in an inverted pattern. However, a perusal of the literature 
shows that this interpretation lacks rigorous proof; its basic assumption - that between 
the two lines observed always the HF line is out-of-phase and the LF line is in-phase - 
has never, to our knowledge, been proved. This stimulated us to take a closer look at 
the whole problem. 

FMR spectra are commonly interpreted on the basis of macroscopic theories; in 
Artman-Layadi theory [IHI^CH]; often applied to bilayer film FMR, sublayer magneti- 
zations are regarded as classical vectors and appear as such in the key expression for free 
energy of interactions between sublayer magnetizations. A microscopic theory of FMR 
in thin films has been developed by some authors, including Ferchmin ^Hl, Puszkarski 
^3 GHj and Wojtczak ^niED]- The Hamiltonian considered in their studies includes 
the exchange interaction energy, Zeeman energy and the uni-axial anisotropy energy, 
and is diagonalized by means of Tyablikov-Bogolyubov method 21 . The resulting 
eigenvalues correspond to spin-wave energies, indicating the resonance line positions, 
and the intensities of these resonance lines can be found from the corresponding eigen- 
vectors. Thus, a full theoretical image of FMR spectrum is obtained by using this 
method. 

In this study, the FMR spectrum in a magnetic bilayer film is investigated as a 
function of the interface coupling and the interface anisotropy, on the basis of the 
theory developed by Puszkarski in ^3 QB] • The Hamiltonian of non-homogeneous thin 
film model is specified in Section El and diagonalized in Section IHl The Hamiltonian of 
ferromagnetic bilayer film is diagonalized in Section HJ Section El refers to the bilayer 
FMR spectra composed of two resonance lines, often reported in experimental studies; 
as we mentioned above, in the most common interpretation the low-intensity line is 
related to an optic mode, and the high-intensity line to an acoustic mode. However, 
the FMR spectra computed in our model prove this interpretation is not always right. 
In Section El we show that symmetric bilayer FMR spectrum does not depend on the 
interface coupling value, and the effect of the interface anisotropy on the critical angle 
appearance and value is studied in Section [71 

2 General planar thin film model 

2.1 Assumptions 

In the planar thin film model, referred to as Valenta model [221 122] i the thin 
film specimen is assumed to be infinite and homogeneous in the directions parallel to 
its surface, but finite and generally inhomogeneous along the surface normal. In the 
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Figure 1: The planar thin film model. 



latter direction, the inhomogeneous thin film can be divided into a number of lattice 
planes parallel to the film surface (see Fig. H]). If all the lattice planes are assumed to 
have identical crystallographic structure (the same plane lattice type and equal lattice 
constant values), all atoms within a plane have identical neighbourhood, and thus are 
in identical physical conditions. A spin is placed in each lattice node. Its position is 
defined by the number of the plane to which the spin belongs (/) and the site vector 
within this plane (j); its nearest neighbours can lie not only in the same plane or 
a neighbouring one, but in farther planes as well (this problem is more thoroughly 
discussed in |JHI2nil2nil2ZI)- The number of plane /' spins being nearest neighbours of 
a plane I spin shall be denoted by zw. Both the spins and the exchange interactions are 
assumed to be equal within a lattice plane (but they can differ from plane to plane). 

Two types of coordinate system shall be used below. In the first one, xyz, related to 
the crystal lattice, x and y axes are parallel to the film surface, the z axis being oriented 
along the surface normal. Another coordinate system to be used, x'y'z', is related to a 
single lattice site; in this local system, the positive direction of z' axis is defined by unit 
vector 7^j, indicating the direction of quantization of the spin in the considered site. 
The direction of quantization follows that of the effective magnetic field in the spin 
site. Three components contribute to this effective field: the applied magnetic field, 
the demagnetization field and the anisotropy field, the latter usually due to the atomic 
magnetic moment interactions with the crystal lattice electric field (the so-called crystal 
field). The assumptions made within this model imply that the effective field within a 
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lattice plane is equal in all spin sites {H^J^-^ = H^^^), and consequently, all the spins 
within a single plane have equal direction of quantization (7^^ = 7z)- 

2.2 The Hamiltonian 

We are going to consider a Heisenberg Hamiltonian of the following form: 

^ = -2 E Jiii'Ai ■ - 9f^B E H!" ■ 4- (%) ' , (1) 

Jijiiji denoting the exchange integral, g being the gyromagnetic ratio, and /i^ being 
Bohr magneton. Vector j specifies the node position in plane / (j G /), and vector 
j' plays the same role in plane /' [j' G I'). Symbol Y^ij-iif means the summation 
involves all the spin pairs, each pair being considered once only. The three terms in 
Hamiltonian (0) correspond to the exchange interactions energy, the Zeeman energy 
and the single-ion anisotropy energy, respectively. The dynamic dipolar fields shall not 
be taken explicitly into consideration here, since - as it has been shown by Krawczyk 
PHj - their effect on the localized state existence conditions is minor. 

Assuming that the exchange interactions occur between the nearest neighbours 
only, and admitting that these interactions, though equal within a lattice plane, can be 
different between spins lying in different planes (in this case, Jij i/j/ = Jw), we obtain 
the following form of the Hamiltonian: 

E Jii'%-^S,j,-gfiBT.H!" ■^S.j-Y.DiiS^j)' ■ (2) 

iJj^i'J' 1,1 ij 

In order to find the inhomogeneous thin film eigenstates, Hamiltonian (j21) must be 
diagonalized. 

3 Hamiltonian diagonalization procedure 

Four transformations shall be performed ^], leading to the diagonal form of Hamil- 
tonian Q: 1) transformation to the local coordinate system, 2) transformation to the 
boson operators, 3) Fourier transformation in the film plane and 4) the transformation 
along the film surface normal. 

3.1 Transformation to the local coordinate system 

First we shall transform Hamiltonian ^ from its form in the crystal lattice coor- 
dinate system, xyz, to that in the local coordinate system, x'y'z' . This transformation 
reads pi] : 

% = ^il^'io + 77f (Aj'S'';; + ^*ijS'i^ , (3) 
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7;j denoting the previously introduced quantization axis unit vector, and Aj^j having 
the following coordinates (provided that y' lies in plane xy): 

h 
h 

Vectors and A^j fulfill the following relations: 

Aj ■ 7ii = 0, a;.. 7,j = 0, 

3.2 Transformation to second quantization operators 

In the second step, the spin operators (expressed in x'y'z' coordinates) shall be 
replaced with boson operators by means of Holstein-Primakoff transformation. In the 
resulting Hamiltonian, only quadratic terms shall be retained, and all the terms of 
other orders shall be omitted (however, the commutation rules should be taken into 
consideration, and attention paid to the sequence of boson operators when omitting 
the higher-order terms). 

Holstein-Primakoff transformation reads pT] : 

= ^ /^jfij^ij^ S' lj = \/2%at/Jj, S'g = Sg - ata^j, 

with operators of creation {ag) and annihilation (dg) satisfying the following commu- 
tation rules: 




J] jj 

Aj ■ Aj 



1, 

0, 



(5) 



Using the approximation of quasi- saturation, we obtain: 

The above transformations shall be applied to the first two terms in the Hamiltonian, 
i.e. the exchange (bi-ion) interaction term and the Zeeman (single-ion, linear) term. 
The transformations to be applied to the remaining anisotropy term (single-ion, non- 
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linear), are as follows ^Uj : 



Ij^ Ij 

Ij^ Ij 



ij"^ ij 
ij'=> ij 

/j*^ Ij 



2% + (25, 



=5;i-(2%-i)25%, 



2%%"«J' 



2%(2%-l)%%, 



2% (25, 



1 atat. 



(9) 



After transformations (jHI), ((HI) and Hamiltonian Q becomes: 
where: 



(10) 



a' i 

~N DiSi (Si cos^ 9 + - cos^ - - 
^ V 2 2 



eff 



TD33 



2 En ZlnJlnSnll " 7n + df^sHf^ ■ Jl 

-Di (5; - I) (1 - 3 cos2 ^) for /J = 



jj' 



Qw 



-D,^Si{Si - |)(1 - cos2 6) for Ij = I'f 
-2y^Ja'ArAi, for Ij ^ I'f 



(11) 



(12) 



(13) 



Rl = -2y^^ ZinJlnSnAi ■ 7„ - giiB\[^iHf^ ■ Ai 

n 

-]-^i{2Si - 1) cos^Vl-cos2^. 



(14) 



Symbol 9, appearing in the above equations, denotes the angle between the film mag- 
netization vector and the surface normal; zw, as previously defined, is the number of a 
plane / spin nearest neighbours in plane V . Equality Jw = Jvi implies that P;-^/ = P^P 
and Qf// = Q^i , which means that Hamiltonian ()1()|1 is Hermitian. 
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3.3 Bi-linear Hamiltonian 

As mentioned above, only quadratic terms shall be retained in the considered 
Hamiltonian: 



The terms of order zero, shifting uniformly the energy scale, are omitted. The first- 
order terms vanish when the direction of 7 corresponds to the minimum energy |HU] . 
The fourth-order terms are related to the interactions between magnons, and their 
omitting is justified by our previous assumptions. 

The general method of quadratic form diagonalization was proposed by Tyablikov 
and Bogolyubov, and applied to thin film and bilayer film by Puszkarski jnBEOl- The 
procedure, analogical to that used by Ferchmin ^| and Corciovei 1^ , is based on two 
Fourier transformations, performed in the film plane and along the direction normal to 
the film surface. By these operations, boson operators and du:f,, originally expressed 

in the direct space, Ij, are transformed into the reciprocal space, k\\k±, in which the 
Hamiltonian is diagonal. 

3.4 Fourier transformation in the film plane 

The transformation in the film plane reads QEl ■ 



denoting the total number of spins in a single lattice plane parallel to the surface, 
and fc|| = [kx, ky] being a vector from the two-dimensional Brillouin zone (coordinates 
kx and ky are quantized through imposing Born-Karman cyclic boundary conditions in 




(15) 




(16) 




(17) 



By transformation ()16|1 Hamiltonian (fTH|) becomes: 



n 



fell,n' 




-hi' 



) 




where: 




(19) 



2 ZlnJlnSnll " In + gj^BHi'''^ ■ - Di (Sl - - j (l - 3 COS^ 6^) 



n 
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(20) 



Term F^^'J, referred to as structural sum, is defined as the following sum over spin 
neighbours: 



and satisfies the following relations: 

^ w — ^ IV ' in — w ^ w — ^11' 



(21) 



(22) 



3.5 Transformation along the surface normal 



In the last step of the diagonalization procedure, canonical Tyablikov-Bogolyubov 
transformation |2II (along the film surface normal) shall be applied to Hamiltonian 



E 



m«(^-l)| 



knk I 



ku,-k I 



(23) 



E 



■kn,-k I 



where and ^j:^^f^^ are operators of creation and annihilation, respectively, of a spin 

wave with energy E{k\\, k±) and wave vector k = [k\\, k±\, components fcy and k± being, 
respectively, parallel and perpendicular to the surface. When expressed by operators 
^■r; , and ft , , the Hamiltonian becomes diagonal: 



n= Eih,k^)^l.L.+ const. 



(24) 



kwM I 



Hamiltonian is, by assumption, Hermitian, which implies that E{k\\, k±) = E*{k\\, kj 
For transformation (j^Hj) to result in the diagonal form of the Hamiltonian, functions 
Ui{k±) and vi{k±) must satisfy the following conditions of orthonormality 



E [Mki.)uiik'^) - vi{k^)v^{k'^)] = 6k^k'^, 

E {Mki.)vi{-k±) - vi{k±)ui{-k±)] = 0, 
I 

E [Mki.)u^'ik±) - vi>{-k±)vl{-k^)] = 5ii' 
E {Mk±K{k^) - v*{-k^)ui,{-k^)] = 0, 



(25) 



and, moreover, they must be solutions of so-called Tyablikov-Bogolyubov equations. 
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In order to find functions ui{k±) and vi{k±) in tlieir explicit forms, we sliall write 
tlie Heisenberg equations of motion for operators 6^ and : 



ib 



kill 



(26) 



fcii fc I 



tk,k^^^ = mhk^)^k,k^^ 



(27) 



By inserting (PHjl into (pUj) . and using (j2Zj), tfie following equation is obtained: 

^ Lik^)Ei^, A;^)^,^ - ^;(_A;^)i5;(_^ll, -A;^)e+, 

kj_ 

EE h'(^ii) U(A:±)4„., + ^^P(-fc±)e!. ) + 

Qr,(fcy) (<'(-A:±)^,,^|^_,^ + t^r(A:±)&„,., 
equivalent to a set of 2L equations, if ^^^^^^and ^ are linearly independent: 

MKfc±)^(^l|, A;^) = E , 

-v!{-k^)E{-k\\,-k^) = J2 [PaihK{-k±) + QUh)uU-kA 



Reversing the direction of k in the second equation 
leads to: 



— fell, ki 



ui{k_i)E{k\\,k±) = J2 [Pii'{h)ui'{ks_) + Q1i,{k\\)vv{kx 
v 

-vi{k^)E{k\\,k^) = [p*,{-k\{)vv{k^) + Qw{-k\{)uv{k, 



-k±) finally 
(28) 



This is the set of Tyablikov-Bogolyubov equations corresponding to our problem. The 
spin wave energy, E{k\\, k±), as well as functions ui{k±) and vi{k±), describing the spin 
wave precession, can be deduced from its solution. 

4 The particular case: ferromagnetic bilayer film 

4.1 General form of the bilayer Hamiltonian 

Let us consider now a magnetic bilayer film, composed of two homogeneous sub- 
layers A and B having identical crystallographic structure; each sublayer is assumed 
to contain a number of lattice planes, Na and A''^, respectively (Fig. 12)). Moreover, 
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Figure 2: The planar magnetic bilayer film model (see the text for detailed description). 

the spins in the lattice nodes, as well as the exchange interactions between the near- 
est neighbours, are assumed to be equal within each sublayer, though they can differ 
between the sublayers. The same rule applies to the other parameters, such as the 
bulk anisotropy constant and the effective field. Besides the standard notion of bulk 
anisotropy, two other quantities shall be used: surface anisotropy D-^^^^ and 

interface anisotropy D^f\ 

The exchange interaction through the interface shall be described by introducing the 
interface exchange integral, Jab, positive for ferromagnetic interface coupling, negative 
for antiferromagnetic interface coupling, and zero in the case where no coupling is 
present between the sublayers. In practice, a thin non-magnetic interlayer is inserted 
between the ferromagnetic thin films A and B, and the interface properties {D^f\ 
Jab) are determined by the interlayer thickness and material. 
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With these assumptions, the matrix form of ()28|) is: 



^AB 

Y 

V 



Y 



\ 



Xa Xab 

^AB J 




2Lx2L 



U and V being defined as follows: 



U 




(29) 









U2 




V2 




, v = 





















(30) 



In films obtained from cubic crystal surface cuts (001), (110) or (111), the nearest 
neighbours of a plane / spin lie in planes up to / + 3 |23 12^1 12Z|; so the general form 
of matrices X and Y is seven-diagonal: 



X, 



X 



B 



Ra — Oo 


Ca 


Da 


Fa 












c\ 


Ra — cii 


Ca 


Da 


Fa 










Da 


C\ 


Ra — 0^2 


Ca 


Da 


Fa 








n 


D\ 


C*A 


Ra 


Ca 


Da 


Fa 










FX 


D*A 


cx 


Ra 


Ca 


Da 


Fa 








FX 


DX 


CX 


RA-h2 


Ca 


Da 










FX 


DX 


CX 


RA-bi 


Ca 












FX 


DX 


CX 


Ra — bo 


Rb — Co 


Cb 


Db 


Fb 












C*B 


Rb - ci 


Cb 


Db 


Fb 










D*B 


C*B 


Rb - C2 


Cb 


Db 


Fb 










D*B 


C*B 


Rb 


Cb 


Db 


Fb 










n 


D*B 


C*B 


Rb 


Cb 


Db 


Fb 








F*B 


D*B 


C*B 


Rb — d2 


Cb 


Db 










F*B 


D*B 


C*B 


Rb — di 


Cb 












Fh 


D*B 


C*B 


Rb — d 



LaxLa 

(31) 



(32) 
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Xab 





Fab 

Dab Fab 

Cab Dab Fab 



LaxLb 



(33) 











^AB 









R\ 





^AB 


^AB 










^A,int 


'^AB 


^AB 


^AB 


{fIbT 


(fIbT 




^B,int 













Rb 









(fIbT 








Rb 



The symbols used above are defined as follows: 



Rb 



^B,s 



LxL 

(34) 





—2SiJiTi, 




A = 


—2SiJiT2, 


(35) 


Fi^ 


—2SiJiTs, 





Ri = -2S,J,To + 9I^bH:" -i. + D.^S, 
+2[z^ + 2{zl + zl + zl)]SiJi, 



R: 



where i — A; B, and: 



,(mt,s) 



-F>i,(int,s)\l SiiSi--] sin^ 



3 cos' 



0i-l) 



(36) 
(37) 





-2^ 


1 SaSbJabA.\ 


AbTi, 




Dab = 




1 SaSbJab-A\ 




(38) 






1 SaSbJab-A\ 


Ab^s, 




^AB — 


-2^ 


1 SaSbJab-^a 


AbT,, 




D\b = 




1 SaSbJab^a 




(39) 


FIb = 


-2\ 


1 SaSbJabAa 


AbTs, 
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ao = 2{z^ + z^ + zf)SAJA'g^iBK^■lA-{Df-DA){SA-\){?>co^''eA-l), 

ai = 2{z^ + zi)SAJA, 

02 = 2z^SaJa, 

62 = 2z^SaJa - 2z^ SbJabIa ■ 1b, 

h = 2{z^ + z^)SaJa-2{z^ + z^)SbJab1a-1b, 

bo = 2{zf + Z^ + Z^)SaJa - 2{Z^ + Z^ + zi)SBJABlA ■ IB - 

gfiBKit ■ lA - {Dit - DaKSa - l){3 cos' 9a - 1), (40) 
Co = 2{z^ + + zi)SBJB - 2{z^ + z^ + z^)SaJab1a ■ IB - 
gfiBKFnt ■ IB - {Dlt - Db){Sb - \){3 cos' 9b -I), 

Cl = 2{Z^ + zi)SBJB-2{z^ + Z^)SAJAB7A-7B, 

C2 = 2zf SbJb - 2z^SaJab1a ■ 1b, 

d2 = 2z^SbJb, 

di = 2{z2 + z^)SbJb, 

do = 2{z^ + z^ + zi)SBJB-9f^BKf -IB- iDf-DB){SB- ^){3cos' 9b- 1). 

As assumed previously, the direction of quantization appearing in the above formulae is 
constant within a single sublayer, but can differ between the sublayers (i.e. 7^ = 7a(b), 
which implies Ai = Aa{b))- 



4.2 The effect of various surface cuts 

We shall henceforth consider a bilayer film whose both sublayers are made of the 
same material {Sa = Sb = S and Ja = Jb = J), possible asymmetry being due only to 
interface (-D^^ 7^ Dfnt^ ^ Kj^^) or surface {K^ ^ Kf) conditions. In the simplest 
case, the nearest neighbours of a plane / spin lie in either the same or a neighbouring 
plane (/' = Z, Z ± 1), which brings the Hamiltonian matrix to a three-diagonal form [23 • 
This situation takes place for surface cut (001) in all three cubic crystal types, and for 
cut (110) in sc and bcc crystals. As our study is focused on ferromagnetic resonance, 
in which standing spin waves only are observed, we assume k\\ = 0. 

In order to simplify the problem, the bilayer spin precession shall be henceforth 
assumed to be circular. The spin precession ellipticity is allowed for in (j29p through 
matrix Y (jHH), whose diagonal and out-of-diagonal elements are defined in (jH7|) and 
(jSni), respectively. If the quantization vector 7 has the same direction in both sublayers, 
then Aa = Ab, and the out-of-diagonal elements of Y vanish, since product Aa ■ Ab 
is zero. The diagonal elements of Y can be divided into three groups, distinguishing 
bulk, surface and interface elements. Their different values are due to the fact that the 
ellipticity of spin precession on the surface and the interface differs from that in the 
bulk. However, if the spin precession ellipticity is assumed approximately homogeneous 
throughout the bilayer film, all the diagonal elements of Y become equal; as shown in 
|30j . this does not affect the relative intensities of the resonance lines, and thus matrix 
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Y, having no effect on the resonance spectrum, can be omitted {i.e. we can assume 
Y = 0). Hence, in the circular spin precession approximation, set of equations 
can be separated, and reduced to a simple eigenvalue problem with matrix X only. 
With the above-specified assumptions, for cubic crystal surface cut (001) equation 
divided by 2SJTi, becomes: 



2 ^^r/ 
-1 


-1 

2 -1 
-1 


2 


-1 










-1 




J int 










Jint 




-1 










-1 


2 -1 

-1 2 -1 

~1 2 — Af^^j- 



U = E'U, 



the Hamiltonian matrix elements being defined as follows: 



(41) 



A 



A{B) 
surf 



1 - 



2SJz 



2SJzi 



(3 cos^ 9 



1 J i.rii 



2SJzi 



0-32 



A(B)^ 
int J 



D 



A{B) 
int 



2SJzi 



(3 cos^ 9 



1x2 



J, 



int 



Jab/ J- 
Properties of function F 



imply that Tui{k\\ = 0) = Z;;/, and thus 
E' = E/{2SJz^), 



(42) 



(43) 



(44) 



(45) 



where zs_ = zi^i+i. 

Considering symmetric surface or interface conditions, we shall henceforth omit 

indices A and B in the respective sublayer parameter symbols, e.g. = = aj2, or 
aA — ab —a 

^surf — ^surf — ^surf- 

Equation ()4ip provides the basis for our further investigation. The Hamiltonian 
matrix corresponds to a simplified bilayer film model on which a spin-wave spectrum 
can be studied as a function of three essential structural magnetic parameters: the 
surface parameter, Ag^rf, the interface parameter, Aint, and the interface coupling, 
Jint- Below, eigenvalue problem ()4ip shall be solved numerically only, assuming the 
spin value equal to one (S=l). 
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5 Ambiguity in existing interpretation of bilayer 
FMR spectra 

In the 'reversed' double-peak FMR spectrum, reported in bilayer films, the high- 
intensity line is commonly interpreted as corresponding to an acoustic mode, the low- 
intensity line being related to an optic mode jHSl IHEl IHH IHHl (HH] • This implies that such 
'reversed' spectrum should appear only in bilayers with antiferromagnetic interface 
coupling. Fig. El shows examples of bilayer resonance spectra with inverse intensity 
arrangement, resulting from our numerical computations; (a)-(c) were obtained assum- 
ing antiferromagnetic interface coupling ^, while (d)-(f) correspond to bilayers in which 
the interface coupling is ferromagnetic. Perpendicular configuration {6 = 0), as well 
as the absence of the uni-directional anisotropy (ojo = 0) were assumed in all cases. 
The case depicted in Fig. El^a) corresponds to the commonly used interpretation: optic 
bulk mode n = 1 has lower intensity than acoustic mode n = 2. In case (b), the 
high-intensity mode (n = 2) is still acoustic, but the low- intensity line (n = 1) corre- 
sponds to an optic interface mode. Another possibility is shown in Fig. Et^c): both first 
modes, n = 1 and n = 2, are of 'optic' nature, n = 2 being an optic bulk mode and 
n = 1 being an optic surface mode; the low intensity of the latter is due not only to its 
optic character, but also to its localization at the surface (this case has already been 
studied in jlO]). If the interface coupling is ferromagnetic, three possibilities can occur, 
their examples shown in Figs. ISl^d)-(f). In case (d), the low-intensity peak (n = 1) 
corresponds to an acoustic surface mode, while the high-intensity line (n = 3) is re- 
lated to an acoustic bulk mode. In (e), the low-intensity mode (n = 1) is acoustic and 
localized at the surface, and the high-intensity mode (n = 2) is acoustic and localized 
at the interface. In (f), the low-intensity mode (n = 1) is acoustic and localized at the 
interface, while the high-intensity peak (n = 3) corresponds to an acoustic bulk mode. 

These results clearly indicate that, from the theoretical point of view, the commonly 
used interpretation of the double-peak resonance spectra, relating the high-intensity 
line to an acoustic mode, and the low-intensity line to an optic mode, is not always 
legitimate. 

6 The effect of interface coupling 

In this Section we shall analyse ferromagnetic resonance spectra obtained in the 
so-called perpendicular configuration, i.e. for 6 = 0. For simplicity reasons, we shall 
assume that there is no uni-directional anisotropy on the interface (ojo = 0), and the 
surface spins have 'natural' freedom {Agurf = 1), Thus, the only source of anisotropy 

^This case (i.e. a bilayer with antiferromagnetic interface coupling) was also considered in our 
previous paper However, as we have just recently realized, the numerical calculations performed 
there - for this particular case only - were incorrect due to the faulty computational program used in 
that paper. This resulted in misinterpretation of the energetically lowest mode as a symmetric one. 
Our present correct numerical calculations show that, in fact, this mode is of antisymmetric nature. 
This discrepancy leads to quantitatively different resonance spectra, however, it does not affect the 
correctness of our main hypothesis formulated in [6'6\ about the permissible interface-localized nature 
of the first resonance mode. 
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Figure 3: Profiles of tlie first {i.e. liaving tlie lowest energies) bilayer modes and the 
corresponding resonance spectra obtained for the parameter values specified in boxes 
(the values of interface parameter A^nt correspond to those assumed for Jj„t and 0^2, 
specified on the right of the bracket). The dashed line indicates the position of a 
hypothetic uniform mode (UM) with A: = 0. 

(a) Asymmetric interface, antiferromagnetic interface coupling. 

(b) Asymmetric interface, interface coupling still antiferromagnetic, but weaker than 
in (a). Note that the value of parameter is negative. 
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Figure 3: c, d 

(c) Asymmetric surfaces, antiferromagnctic interface coupling. 

(d) Asymmetric interface, ferromagnetic interface coupling. 



18 



^surf 



^int 



= 1 

= 0.9 



Aint^ = 0.9 



^surf 



= 1.5 



(e) 

ai2^=-0.1 

Ji„t=0.3 

ai2^=-0.1 



UM 
2 



5 43 
J LX 



ENERGY 




n=3 



n=2 



n=l 



n=l 



^surf 



= 0.5 



Aint — 1.3 



Aint — 1.3 



Asurf = 0.6 



(f) 

ai2^=-0.3 

Ji„t=0.3 

ai2^=-0.3 



7 
I 



UM 



ENERGY 



11=4 




n=2 



Figure 3: e, f 

(c) Asymmetric surfaces, ferromagnetic interface coupling, 
(f) Asymmetric surfaces, ferromagnetic interface coupling. 
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Figure 4: Energies of symmetric bilayer modes versus the interface coupling integral, 
Jint- The solid lines correspond to antisymmetric modes, the dotted lines corresponding 
to symmetric ones. 

is the interface uni- axial amsotTopj. Fig. 0] shows bilayer mode energies as functions 
of the interface coupling integral in the absence of any interface anisotropy (0^2 = 0) 
(this situation was studied in |SS!)- If the couphng is ferromagnetic, Jint has no effect 
on the energy of the odd modes, while in the case of antiferromagnetic coupling, the 
modes 'insensitive' to variations are even. This result fully corresponds to that 
reported in jHS] and [201 • An insight into these functions is provided by the mode 
profiles analysed versus the interface coupling integral. Fig. Efa) shows profiles of the 
six lowest modes for five different Jmt values. In the case of ferromagnetic coupling the 
interface coupling integral is found to have little effect on the shape of these profiles. 
Mode n = 1 is always a uniform mode; the even modes are antisymmetric (and as 
such do not appear in the resonance spectrum), and all the odd modes except n = 1, 
though symmetric in the bilayer, are antisymmetric within each sublayer and thus do 
not appear in the resonance spectrum either. Hence, it is only the uniform mode that 
is observed in the resonance spectrum (see Fig. E^c)). When J^nt < 0, the lowest 
mode is localized at the interface (the localization becoming stronger as Jint absolute 
value increases), while mode n = 2 is uniform. However, also in this case only one 
resonance line is observed (n = 2, see Fig. El^b)), as all the odd modes, including the 
interface mode, are antisymmetric, and thus do not appear in the resonance spectrum. 
Hence, the resonance spectrum is found to be insensitive to the interface coupling 
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Figure 5: (a) Profiles of tlie six lowest modes and (b, c) the corresponding FMR 
spectra in a symmetric bilayer, computed assuming no interface anisotropy (0^2 = 0), 
the surface spins having natural freedom {Ag^rf = 1), for different values of interface 
exchange integral Jint- Each sublayer is composed of = 11 lattice planes. UM 
denotes the resonance line corresponding to the uniform mode excitation. 
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integral variations (since Jint has no effect on the symmetric modes, and affects only 
the antisymmetric ones). 

The above conclusion can be equally deduced in a different reasoning, based on 
the notion of effective interface parameter, A^ff, introduced in [3H|. In the considered 
perpendicular configuration {6 = 0) this parameter is expressed as follows: 

^ _ i Alj f = 1 — 2aj2, for symmetric modes, . , 

'^■^^ I Agjj = 1 — 2aj2 — 2 Jint, for antisymmetric modes. 

The above formulae directly indicate that only the antisymmetric modes depend on 
Jint, the symmetric modes being unrelated to it. This remains valid also when 7^ 0, 
as shown in Figs. IHl and [7| However, the FMR spectrum depicted there contains more 
than one peak, the resonance lines corresponding exclusively to either odd or even 
modes, for Jint > or Jint < 0, respectively. In each previously, the FMR 

spectrum does not depend on Jint value. 



7 Critical angle in bilayer resonance spectrum 

The resonance spectra considered in the preceding sections were obtained for the 
so-called perpendicular configuration, in which the magnetization vector is oriented 
along the film surface normal. In this section, we are going to investigate configuration 
effects in resonance spectrum, due to different orientation of the magnetization vector 
with respect to the film surface. Our analysis shall be based on the effective interface 
parameter, Agff, a concept introduced in our earlier study jHHl- In a symmetric bilayer 
film, this parameter is expressed as follows: 




Alj-j = 1 — fljo ~ ^42(3 cos^ 9 — 1) , for symmetric modes, 

A\. = 1 — ttjo — cti2(3 cos^ 9 — 1) — 2 Jint, for antisymmetric modes. 



(47) 

Bilayer resonance spectra are generally composed of several resonance lines (so- 
called spin wave resonance - SWR) , though in certain conditions the spectrum reduces 
to a single peak. This occurs, for example, when angle 9 between the magnetization 
vector and the surface normal takes a particular value, referred to as critical angle. The 
existence of the critical angle is due to the fact that the effective interface parameter 
fl47|) is a function of 9. For simplicity reasons, in our investigation of the critical angle 
effect we shall assume that the surface parameter, Asurf, does not depend on 9, and 
that Asurf = 1 (^-6. the spins have 'natural' freedom on both surfaces). 



7.1 The effect of uni-axial anisotropy 

In the first case to be considered, the effect of the uni-directional anisotropy shall 
be neglected (ojo = 0). Figs. |H^a) and|Hfc) show the intensities of mode excitations 
(in relation to that of the first symmetric mode, n = 1) as functions of angle 9 in 
a symmetric bilayer film (composed of 22 atomic planes, 11 planes in each sublayer) 
with ferromagnetic interface coupling {J^nt = 0.3). Graphs (a) and (c) were plotted 
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Figure 6: (a) Profiles of the six lowest modes and (b, c) the corresponding FMR spectra 
in a symmetric bilayer, computed assuming interface anisotropy 0^2 = 0.2, the surface 
spins having natural freedom {As^rf = 1), for different values of interface exchange 
integral Jint- Each sublayer is composed of = 11 lattice planes. BM denotes a 
resonance line corresponding to the symmetric bulk mode excitation. 
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Figure 7: (a) Profiles of the six lowest modes and (b, c) the corresponding FMR 
spectra in a symmetric bilayer, computed assuming interface anisotropy = —0.1, 
the surface spins having natural freedom {Ag^rf = 1); for different values of interface 
exchange integral J^nt- Each sublayer is composed of = 11 lattice planes. IM denotes 
the resonance line corresponding to the symmetric interface mode excitation. 
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Figure 8: Resonance line relative intensities versus the configuration angle in a bilayer 
film with symmetric interface boundary conditions and ferromagnetic coupling {Jint = 
0.3). The other parameter values assumed: = 11, ajo = 0, Agurf = 1- Graphs (a) 
and (c) show intensity ratio In/h plotted versus angle 6; the corresponding profiles of 
the six lowest modes are shown in (b) and (d) for three different 6 values. 
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assuming = 0,^2 = —0.3, and 0,^ = 0.^ = 0.3, respectively. The corresponding 
profiles of the six lowest modes are depicted in Figs. Efb) andlHfd), for three different 
values of angle 9: 40" {9 < 9crit), 54.7" {9 = 9crit) and 70" {9 > 9crit)- (Note that 
in both cases the even modes, being antisymmetric, do not appear in the resonance 
spectrum) . 

As shown in Fig. IH^b), for ~ '^^ ~ —0.3, mode n = 1, localized at the interface 
when 9 < 9crit, becomes a uniform mode for 9 = 9crit, and a bulk mode when 9 > 9crit- 
The other modes are of bulk nature in all 9 range. The resonance intensity of the first 
bulk mode (n = 3), surpassing that of the interface mode at 6^ = 0° (see Fig. lEfa)), 
decreases with growing 9 and equals the interface mode intensity at ^ ^ 23.4°. As 9 
continues to increase, the intensities of the bulk modes decrease still further to vanish 
completely at = 9crit{^ 54.7°), when mode n = 1 becomes uniform, and the only one 
to appear in the SWR spectrum. When 9 > 9crit, all the modes are of bulk character, 
the intensities of modes n > 3 being much lower than that of mode n = 1. Thus, the 
critical angle is found to separate two regions in which the relative mode intensities 
are radically different. 

For = = 0.3 the situation is reversed (Figs. IHl^c) and|Hl^d)). When 9 < 9crit, 
all the modes are of bulk character; mode n = 1 becomes uniform at ^ = 9crit (the 
intensities of all other modes being zero), and localized at the interface when 9 > 9crit, 
its intensity remaining the highest in all 9 range. This case is a 'mirror image' of that 
depicted in Figs. ID^a) andlH^b); also here, the critical angle separates two regions in 
which the relative mode intensities are completely different. 

When the interface coupling becomes antiferromagnetic {Jint = —0.3, see Fig. 
the angle relations of the relative mode intensities (with respect to the intensity of the 
first symmetric mode, n = 2) are in principle similar to those obtained in the case of 
ferromagnetic coupling, the critical angle remaining 54.7°. 

From the numerical analysis presented above we deduce that in the case of symmet- 
ric bilayer with surface spins having natural freedom, the critical angle value is 54.7° 
and does not depend on either the interface coupling, Jint, or the uni-axial anisotropy, 
ai2- Its independence of the interface coupling is a consequence of the fact that Jint, 
having no effect on the symmetric modes, does not affect the bilayer resonance spec- 
trum, as shown in Section IHl As regards the uni-axial anisotropy, its variations, though 
significantly modifying the SWR spectrum, do not change the 9crit value, as the critical 
angle condition, Aljj- = 1, is satisfied when cos^ — 1 = (or 6^ = 54.7°) (jT7|) . being 
uninvolved. 

Let us now consider a bilayer film with asymmetric interface conditions. Fig. [TUl 
shows the corresponding relative mode intensities plotted versus the configuration an- 
gle, 9. In the considered asymmetric bilayer film, each sublayer is composed of 11 
atomic planes, the interface coupling is antiferromagnetic {Jint = ~0.3), and the as- 
sumed interface uni-axial anisotropy values in sublayers A and B are = —0.3 and 

= or 0.2, respectively. 

Asymmetric interface conditions are found to have a significant effect on the res- 
onance spectrum. As the asymmetry becomes stronger, the relative intensities of the 
quasi-antisymmetric modes increase, while those of the quasi-symmetric modes de- 
crease, but the critical angle value remains unchanged and equal to 54.7°, as in the 
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Figure 9: Resonance line relative intensities versus the configuration angle in a bi- 
layer film with symmetric interface boundary conditions and antiferromagnetic cou- 
pling {Jint = —0.3). The other parameter values assumed: = 11, ajo = 0, Ag^rf = 1- 
Graphs (a) and (c) show intensity ratio In/h plotted versus angle 9; the corresponding 
profiles of the six lowest modes are shown in (b) and (d) for three different 6 values. 
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Figure 10: Resonance line relative intensities versus the configuration angle in a bilayer 
film with asymmetric interface boundary conditions and antiferromagnetic coupling 
(a^ 7^ a^, Jj„t = —0.3). The other parameter values assumed: = 11, ajo = 0, 
Asurf = 1- Graphs (a) and (c) show intensity ratio In/h plotted versus angle 9 for two 
different pairs a^, a^; the corresponding profiles of the six lowest modes are shown in 
(b) and (d) for three different 6 values. 
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case of symmetric bilayer film. This is due to the fact that in the critical configuration 
a bilayer film must be fully symmetric. 

Thus, in a ferromagnetic bilayer film with 'natural' surfaces and the interface prop- 
erties defined by two parameters only, namely the interface coupling and the interface 
uni-axial anisotropy, the critical angle value is always 54.7°, for both symmetric and 
asymmetric interface conditions. 

7.2 The effect of uni-directional anisotropy 

The effect of the uni-directional anisotropy on the critical angle value in symmetric 
bilayer SWR spectrum is illustrated in Fig. ^2 Each sublayer contains 11 atomic 
planes, and the coupling between them is of ferromagnetic nature {Jint = 0.3). The 
assumed values of interface uni-axial and uni-directional anisotropics are = —0.3 
and Ojo = 0.2; —0.2, respectively. 

The critical angle is found to strongly depend on the uni-directional anisotropy. 
The exact relation, deduced from (jTf|) . is as follows: 



The above relation is plotted in Fig. (bold line); in the same graph, intensity ratio 
/3//1 is plotted versus angle 9 and uni-directional anisotropy Ojo- Function follows 
exactly the line along which the value of function /3//i(^,aio) vanishes. 

Thus, a critical angle value different from 54.7°, found experimentally in a sym- 
metric bilayer SWR spectrum, would provide evidence for co-existence of uni-axial and 
uni-directional anisotropics in the studied sample, the ratio of these two anisotropy 
types being deducible from the critical angle value (according to (@HI))- 

8 Conclusions 

The results of our theoretical investigation of bilayer FMR spectrum, presented in 
this study, deny the interpretation commonly used in experimental studies reporting 
double-peak FMR spectra in bilayer films. In this interpretation the low-intensity peak 
is related to an optic mode, and the high-intensity line to an acoustic mode. However, 
our theoretical study shows this is not always true, providing examples of double-peak 
spectra in which both appearing modes are acoustic or, what is more, the high-intensity 
peak corresponds to an optic mode. Thus, the interpretation of the observed modes 
cannot be unambiguous a priori, as the intensity of a resonance line depends not only 
on the phase shift in the sublayer magnetization precession, but also on the precession 
amplitude distribution, and especially on its localization. 

In the investigated symmetric bilayer SWR spectra both the position and the in- 
tensity of the resonance lines are found to be independent of the interface coupling 
integral value. The reason is that in symmetric bilayer film the effective interface pa- 
rameter depends on the interface coupling integral for antisymmetric modes only, and 
these modes do not appear in the SWR spectrum. Hence, symmetric bilayer film seems 
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Figure 11: Resonance line relative intensities versus the configuration angle in a bilayer 
film with symmetric interface boundary conditions and ferromagnetic coupling, assum- 
ing non-zero uni-directional anisotropy [Jint = 0.3, aio 7^ 0). The other parameter 
values assumed: N = 11, Agurf = 1- Graphs (a) and (c) show intensity ratio /„//i 
plotted versus angle 9 for two different values of ajo; the corresponding profiles of the 
six lowest modes are shown in (b) and (d) for three different 6 values. 
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Figure 12: Intensity ratio /3//1 of the two lowest symmetric modes (n = 1, 3) in a sym- 
metric bilayer SWR spectrum, plotted versus the configuration angle (angle between 
the film magnetization and the surface normal), 6, and the interface uni-directional 
anisotropy, ajo- The ferromagnetic interface exchange coupling is assumed, with inter- 
face uni-axial anisotropy value 0^2 = —0.3. The bold line represents the critical angle, 
Ocrit, as a function of the interface uni-directional anisotropy, Cjo- 
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to be particularly convenient for interface anisotropy studies. However, it should be 
remembered that when the interface coupling is ferromagnetic the odd modes are sym- 
metric, and the even modes are antisymmetric, while in the case of antiferromagnetic 
coupling the situation is reversed: the odd modes (including the first one, with the 
lowest energy) are antisymmetric, the even modes being symmetric. 

Moreover, we show that the critical angle effect can occur in a bilayer film even 
in the absence of surface anisotropy. In this case, this effect is totally due to the 
interface, and if the interface anisotropy is purely uni-axial, the critical angle value is 
always 54.7° and does not depend on the anisotropy value. Any deviation of the critical 
angle value from 54.7° involves an additional source of interface anisotropy, namely the 
uni- directional anisotropy. 
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Appendix: Bilayer Hamiltonian matrix elements 

Below we specify the exphcit form of Hamiltonian matrix elements for cubic crystal 
surface cuts (001) and (110). Our formulae were derived for k\\ = 0, with the quanti- 
zation vectors in both sublayers assumed to be identical {i.e. 7^ = 7b, which implies 
A A = Ab), and the spin precession to be circular. On these assumptions matrices (jHH), 
and dSSl) become: 



a C A 



Ra C / 



Ca Ra Ca 
Ca RA — b 



X, 



Rb — c Cb 

Rb Cb 



C*B 



Cb Rb 



C 



B 



Cb 
Rb — d 



X 



AB 



... 

Cint ... 



In all the relations detailed below index i denotes the sublayer label (A or B) 
sc(OOl) surface cut 



Ri 

Cint 



4S,J, + gfip{H^^^ ■^) + D,[^S,--]{3 cos' 9 - 1 



2\J SASBJint) 



1, 



2SaJa - QfiBKt ■ 7 - (Df - Da){Sa - ^)(3 cos^ 9-1), 
2SaJa - 2SBJ^nt " g^^BKt " 7 " (A^ " Da){Sa - ^)(3 cos^ 9 - 1) 
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C = 2SBJB-2SAJ^nt-g^iBKfnt■1-iDZt-DB){SB--){3cOs''0-l), 

d = 2SBJB-gi^BK^-l-{Df-DB)iSB-l){^cos''e-l). (49) 
bcc(OOl) surface cut 

Ci — —SSiJi, 

Cint — ~^\/ SASBJint, 

a = SSaJa- QI^bK^ -7 -iDf-DA){SA -^){3cosH-l), 

b = SSAJA-SSBJint-9f^BKf^t-l-{Dtt-DA){SA-l){^cos'9-l), 

C = 8SBJB-8SAJint-gi^BKlt-1-iDfnt-DB){SB-l){3cOS^e-l), 

d = SSbJb- 9I^BK^-l-(Df-DB)(SB- l){Scos''e-l). (50) 
fcc(OOl) surface cut 

Ci = SS^Ji, 

Cint = — SASBJint, 

a = 85^Ja-5)"b^^7-(^.''-^a)(^a-^)(3cos2^-1), 

6 = 8>5^JA-8>5BJ,„t-^/iBi?i^,-7-(At-^A)(5A-^)(3cos'^-l), 

c = 8>SBJB-85^J,„t-^/iB^^,-7-(A^,-L'B)(5B-^)(3cos2^-l), 

d = SSbJb- 9I^bK^ -7- {Df -Db){Sb-^){3 cos^e~l). (51) 

sc(llO) surface cut 

Ri = 85iJ, + ^/.;3(^r^^-7) + A(5.-^)(3cos'e-l), 

Ci = —ASiJi, 

Cint = ~^\J SASBJinti 

a = 4,S^J^-^/iB^,^.7-(L>,^-L>^)(,5^-i)(3cos2^-l), 

h = ASAJA-^SBJint- gi^BKit-l -{Dit- DA){SA-\){2,COs'e -l\ 
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C = 4SBJB-^SAJ^nt-g^iBKfnt■1-iDZt-DB){SB--){3cOs''0-l), 

d = ASBJB-gi^BK^-l-{Df-DB)iSB-l){^cos''e-l). (52) 
bcc(llO) surface cut 

Ri = ASiJ, + gf,p{Hff ■^) + Di (^S,-^^ {3 cos-" 0-1), 

Ci — —ASiJi, 

a = ASAJA-gi^BK^-l-{Dj-DA){SA-\){2>cos^e-l), 

b = ASAJA-^SBJint-gpiBKL-l-(Dit-DA){SA-l){ScOs'9-l\ 

C = 4SBJB-^SAJi„t-gfiBK^nfl-{D^nt-DB){SB-\){^COs'e-l), 

d = ASB-JB-giiBKf -^-{D^ -DB){SB-\){3cos'e-l). (53) 



36 



